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Strong-coupling expansions, to order (t/J)^, are derived for the Kondo lattice model of strongly 
correlated electrons, in 1-, 2- and 3- dimensions at arbitrary temperature. Results are presented for 
the specific heat, and spin and charge susceptibilities. 



I. INTRODUCTION 



This paper, the second of a sequence, studies the thermodynamic properties of the Kondo lattice model, described 
by the Hamiltonian 



The first term describes a single band of conduction electrons, the "Kondo coupling" term represents an exchange 



conduction electron density via a chemical potential. 

The Kondo lattice model combines two competing physical effects. In the strong coupling (large | J|) limit, the 
conduction electrons will form local singlets (J > 0) or triplets (J < 0) with the localized spin at each site. In 
either case there will be a gap to spin excitations and spin correlations will be short ranged. On the other hand, at 
weak coupling, the conduction electrons can induce the usual RKKY interaction between localized spins, leading to 
magnetic order. 

Despite the apparent simplicity of the model, no exact results are known, either at T = or at finite temperatures, 
in any spatial dimension. In the preceding papeiE we studied the ground state properties of the model at T = 0, 
using linked-cluster series expansions. We refer to that paper for a discussion of other work, which has been, almost 
exclusively, restricted also to T = 0. In the present paper we focus on finite temperature thermodynamic properties. 
We know of only a few previous studies of this kind. Roder et a/.i considered the ferromagnetic model in the limit 
J —>■ — oo, on the simple cubic lattice, via a high temperature expansion. Here we treat the general J case and focus 
on the antiferromagnetic model. Shibata et aM have studied the ID antiferromagnetic model via a finite-temperature 
DMRG approach. Haule et aL0 have treated the 2D case, primarily via a numerical finite temperature Lanczos 
method. We compare our results with this work wherever possible. Haule et al. have also considered the atomic limit 
{t — 0), and the order t'^ correction terms. Our work was largely motivated by this paper. 

Our approach, which will be described in the following Section, treats the single site terms exactly and treats the 
hopping term perturbatively. It is, thus, an expansion about the "atomic limit" . We summarize here, for completeness 
and for later reference, the exact results in this limit. 

For variable conduction electron density there are 8 states per site: 2 states with no conduction electrons and 
localized spin up or down, 2 states with two conduction electrons of opposite spin, and 4 states (a singlet and 3 
triplets) with one conduction electron coupled to the localized spin. For a lattice of N sites the grand partition 
function is 




(1) 



interaction between conduction electrons and a set of localized S = ^ spins, and the final term allows for variable 



Zo = Tr{e-''«} 




(2) 
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with 

r,3K I Q„-A' 



zo = 2 + (e^" +3e-")C + 2C' (3) 



and K = pj/4, ( = e^^, /3 = l/ksT. 
The internal energy per site is given by 



Inzo 

C 



2+(e3^ + 3e-^)C + 2C2 
The fugacity can, as usual, be eliminated in favour of the electron density n by using the relation 

d 



Inzo 

T 

where we have introduced q = ■^{e^^ + 3e~^). Solving this gives 



(5) 



-n) + ^q^l - ny + 4n(2^ 



The specific heat is then obtained from (^) and via the usual relation C„ = 
Also of interest are the compressibility or "charge susceptibility" 

dn 

which can be expressed as 



(7) 



and the magnetic susceptibility, which is given by 



1 + gC + C 



1 + + c 

Simpler analytic results can be obtained in various limits: at high or low temperature, and at or near half-filling. 
Some of these are given in Ref. ^, although they contain minor errors in a few cases. 

Of course in the atomic limit all of these quantities are smooth functions of temperature and electron density. 



II. THERMODYNAMIC PERTURBATION THEORY 

Our goal is to obtain an expression for the thermodynamic potential, and hence other quantities, in powers of </ J. 
We work in the grand canonical ensemble and write the Hamiltonian as 

H = Ho + V (10) 

with 

Ho = Si • - nicr (11) 

i ia 
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The grand partition function can then be expanded in the usual way as 

Zg = Tr{e-'^(^«+^)} 

^ Zoll+y^i-iy dn--- dTr{V{Ti)---V{Tr))\ (13) 

r=l "'0 "'0 



with 



and 



V{t) = e^^^Ve-^"" (14) 



(A) = ^Tr{e-''^M} (15) 

^0 



where Zq is the atomic hmit partition function (Q). The free energy (grand potential) is then given by 



PF ^ Nlnzo + Y^i-iyTr (16) 



where 



- dn • • • / " ' dTr{V{Ti) ■ ■ ■ V{Tr))N (17) 

and the subscript N signifies that only the part proportional to N is to be included. 

This approach is, of course, well known and has been used in the past for both pure spin modelsS and for the 
Hubbard model!. Any contribution to Tr in ([l^ comes from a particular cluster of sites and bonds, a "graph" . It 
is possible to restrict the class of graphs to connected ones only, as done in Ref. |[ In our work we included also 
disconnected graphs, of more than one component, as they are rather easy to deal with directly. Since each bond 
contributes a V operator and hence a factor t, it is obvious that to carry the expansion to order all topologically 
distinct graphs with up to r bonds need to be considered. There are a total of 115 graphs through 8th order, which 
is as far as we have been able to compute. 

The contribution of a particular graph G to the free energy can be expressed in the form 

TAG) = Ccit/JTz^P J2 «s,L™i^'e™^C (18) 

s,l,m 

where Cg is the embedding factor, or "weak lattice constant" of graph G in the particular lattice considered, p is the 
number of points or vertices in the graph, the as^i^m are numerical constants, and the sum contains, for each graph, a 
finite set of terms labeled by integers s, /, m. Evaluation of these expressions is a lengthy procedure, involving a trace 
over a space of 8^ states and evaluation, for each term in the trace, of an r-fold multiple integral. It is possible to 
find many time saving refinements, but, even so, the evaluation of the worst case, the octagon, took some 260 hours 
of CPU time on a IGHz Compaq alpha processor. 

Having computed the T(G) factors for all graphs it is then a simple matter to combine these and to obtain, for any 
lattice, the free energy per site in the form 

OO 

-Pf = Inzo + z^^FriK, CKt/jy (19) 

where the Fr are complete expressions in K and C, of the form 

F,(if,C)= 5] as,i,™i^'e™V (20) 

s.l.m 
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These expressions are too lengthy to display here, but can be supplied on request. To give some idea of the size, the 
8th order factor Fg contains 1042 separate terms. 

From ( p^ and ( po|) one can compute expressions for the internal energy and specific heat. The internal energy can 
be expressed in the form 

OC 

= uo + Y,^o^^^'^MK,CKt/jy (21) 

where uq is the atomic limit result (|^) and the Er{K,Q are complete expressions. The specific heat is Cy = du/dT. 
For most purposes it is more useful to express the series in terms of electron density n, which can be obtained from 
(|l|) via 

n = C^(-/3/) 



= no + Y,^o''''^'^YriKX)it/JY (22) 



r=2 

where no is the atomic limit result and the Yr{K, C) are, again, complete expressions in K, C,. For fixed n and K 
we then use a numerical reversion procedure to obtain expansions for the fugacity, in power oit/ J , which can then be 
eliminated from the thermodynamic functions. For half- filling (n = 1) this reversion can be carried out analytically. 
In addition we have included a magnetic field term 

H' = H- hY,{n^,t - n,^i + 25f ) (23) 
to allow calculation of the zero-field magnetic susceptibility Xs, which is expressed in the form 

oc 

p-'Xs = r'xo + E ^o^'^^'^XriK, CKt/jy (24) 

r=2 



where, again, the first term is the atomic limit result (j^) and the Xr are complete expressions. The fugacity can again 
be eliminated in favour of the electron density, as described above. 

Several checks on the correctness of our results have been made. Since temperature enters explicitly, we can take 
the zero temperature {K — ^ oo) limit analytically, to recover the ground state energy serieJiB. Another test is the 
J limit, in which Fr{K, Q in ( p^ has a leading term of order K^, giving an expression in powers of (/3t). This can 
be compared with the results for free electrons. Complete agreement is obtained. We are confident that our results 
are correct, but cannot exclude the possibility of errors not picked up by these checks. 

The analysis of the results follows standard lines. For any fixed K,( we obtain expressions for thermodynamic 
quantities as a series in the single variable {t/J). For small t/J the series converges rapidly and even a naive sum 
gives good accuracy. Fade approximants and integrated differential approximants§ allow extrapolation to larger (t/J). 
For fixed n we obtain series for ( in terms of t/J, which arc then substituted into corresponding expressions (|2C|), 
(^) to obtain a single- variable series in {t/J), which is then evaluated as before. To calculate the specific heat it is 
necessary to include derivatives of ^ with respect to temperature. For the benefit of the reader we give in Table I 
coefficients for some representative cases. Other series can be supplied on request. 

While the general expressions are too lengthy to write down, at high temperatures, one can expand the various 
quantities in powers of (3 or K. The expressions up to order K'^ are provided in the Appendix. These are valid for all 
loosed-packed/bipartite lattices. 
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FIG. 1: The specific heat Cv (a), magnetic susceptibility Xs (b) and charge susceptibility Xc (c) versus T/J for the linear chain 
at n = 1 for t/J = 0, 0.2, 0.4, 0.5, 0.6. 



In the following sections we present results for the linear chain, the square lattice (sq) and the simple cubic (sc) 
lattice. A discussion of other lattices will be presented elsewhere. 



III. THE ID KONDO LATTICE MODEL 

We consider the antiferromagnetic model on a linear chain, with arbitrary electron concentration n. Using the 
procedure described above we have computed the specific heat and spin and charge susceptibilities. Figure 1(a), 
(b) (c) shows these quantities, as functions of temperature, for n = 1 (half-filling) for t/J = 0,0.2,0.4,0.5,0.6. For 
t/J < 0.4 the series are well converged and the curves are obtained simply from the partial sums. For the larger values 
integrated differential approximants have been used, with the error bars indicating the variation between different 
approximants. Good agreement is obtained with the DMRG results for all three quantities. However our (relatively 
short) series are unable to probe the larger t/J region as well as DMRG can. 

Some comments on Fig. 1 are in order. Fig. 1(a) shows the specific heat, which shows an interesting crossover 
from a single peak for small t to a two-peak structure at larger t. In the larger t region the high temperature peak 
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becomes broadened and less prominent and is, presumably, due to conduction electrons whereas the low-temperature 
peak arises from the fluctuating local spins. The spin susceptibility (Fig. 1(b)) also has a peak at a characteristic 
temperature. The peak is enhanced and moves to lower T on increasing the hopping parameter t. Increasing t/ J will 
weaken the singlet correlations and a lower characteristic temperature is sufficient for thermal fluctuations to become 
dominant. There is some indication of a double peak for t/J = 0.6, but this may be an artifact of the numerical 
analysis. The charge susccptibiUty also has peak at a characteristic temperature, and a rapid drop to zero at low 
temperatures. The peak is depressed with increasing t but the position stays relatively constant. 

To display the effect of varying conduction electron density we have chosen an intermediate value t/ J = 0.4 and 
show, in Fig. 2(a), (b), (c) curves of C„, Xs and Xc versus temperature for n = 1, 0.75, 0.5 and 0.25. The high- 
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FIG. 4: The specific heat C„ versus T/J for square lattice (a) and simple cubic lattice (b) at t/J = 0.3 for n = 1, 0.75, 0.5, 0.25. 
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FIG. 5: The magnetic susceptibility X" versus T/J for square lattice (a) and simple cubic lattice (b) at n = 1 for t/J 



temperature peak in Cy (Fig. 2(a)) drops roughly proportionally to n, in agreement with the assignment of this peak 
to conduction electrons. The series do not allow the low-T specific heat to be determined with sufficient precision to 
see the effects of doping. Tsunetsugu et aM present a ground state phase diagram of the ID Kondo lattice (Fig. 6 of 
Ref. where a transition line separates a small n ferromagnetic phase from a large n paramagnetic phase. For the 
parameter ratio t/J = 0.4 the critical doping is ric ~ 0.65. While there cannot be true order at finite temperature, 
the marked change in the low T specific heat between the two curves with n = 0.75 and 0.5 may will be a reflection of 
this effect. The effect of doping on Xs is dramatic. For n = 1, at low temperatures, the system is in a gapped singlet 
phase and Xs goes to zero exponentially. Away from half-filling there will be free spins and Xs diverges according to 
the usual Curie law. The charge susceptibility Xc also shows a sharp crossover between the undoped and doped cases. 

The ID results are consistent with the, presumably, more accurate DMRG calculations. We are not aware of any 
published DMRG results for specific heat at finite doping, or for the susceptibility for n < 0.8. Our results confirm 
that the series approach can be successfully applied to this model. We now turn to the 2D and 3D cases, where far 
less is known and where other methods have particular difficulties. 



IV. THE 2D AND 3D SYSTEMS 

We have computed and analysed series for the specific heat, spin susceptibility and charge susceptibility for both 
the square and simple cubic lattices. Some representative series are given in Table I. We choose to present results for 
the two lattice together so as to highlight similarities and differences between them. 

Figure 3 shows the specific heat at half-filling for various t/J ratios. The qualitative behaviour is similar to the 
ID case, although smaller values of t/J suffice to produce comparable deviations from the atomic limit. A two-peak 
structure is manifest a±t/ J — 0.4. Figure 4 shows the dependence of the specific heat on doping (for fixed t/ J — 0.3). 
The decrease with decreasing n again confirms that the high T specific heat is due to conduction electrons. Figures 5 
and 6 show the magnetic susceptibility, at half-filling for various t/J and at t/J — 0.3 for various electron densities. 
Figure 6 shows a stricking crossover from n = 1, where Xs is structureless and vanishes as T — » to lower densities 
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n = 0.5 and 0.25 where Xs appears to diverge. Evidently n = 0.75 is near the critical concentration, as Xs appears to 
start to diverge but then drops to zero. There is also some apparent structure for the 3D case at T/J ~ 0.15, which 
may be an artifact of the analysis. 

A mean-field treatment! suggests the existence of a finite-temperature phase transition to a ferromagnetic phase 
for small n, at least for the 3D case. To explore this we plot the inverse susceptibility versus T in Figure 7. For both 
lattices x7^ appears to vanish linearly as T — > at n = 0.5. For the simple cubic lattice there is some indication that 
the curves for xj^ aX n = 0.25 and n = 0.1 vanish at a small finite T, confirming a finite temperature ferromagnetic 
phase. However the analysis at low temperature is imprecise and does not allow for an accurate determination of Tc. 
The data for the square lattice appear to show no finite temperature ferromagnetic transition. 

Finally in Figures 8 and 9 we show the charge susceptibility versus temperature. The behaviour is qualitatively 
similar to the ID case. 



V. CONCLUSIONS 

We have used thermodynamic perturbation theory to investigate the antiferromagnetic Kondo lattice model at 
finite temperatures. Our calculations focus, in particular, on the specific heat and spin and charge susceptibilities, 
and their variation both with the ratio t/ J and with electron concentration. We have presented results for the linear 
chain, square lattice, and simple cubic lattice. 

Overall, for the parameter region where our series can be analysed with reasonable precision, the behaviour of the 
3 lattices is qualitatively similar. This is to be expected at moderate and high temperatures. We do, however, see an 
indication of a finite temperature ferromagnetic transition in 3-dimensions for small n, consistent with expectations. 
This is not seen in 1- and 2-dimensions. 

For the linear chain our results are in excellent agreement with the previous finite-temperature DMRG calculationsi, 
and serve to confirm the accuracy of both methods. There are few existing results in higher dimension and our work 
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FIG. 6; The magnetic susceptibility Xs versus T/J for square lattice (a) and simple cubic lattice (b) at t/J = 0.3 for n 
1,0.75,0.5,0.25. 
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FIG. 7: The inverse of magnetic susceptibility 1/xs versus T/J for square lattice (a) and simple cubic lattice (b) at t/J = 0.3 
forn = 0.75, 0.5, 0.25, 0.1. 



should provide a valuable benchmark for other approaches. While we have not attempted a detailed comparison 
on or fit to experiment, this would be possible. We plan to report results for closed-packed lattices (triangular and 
face-centered cubic) and for the ferromagnetic Kondo lattice model elsewhere. 
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FIG. 8: The charge susceptibility Xc versus T/J for square lattice (a) and simple cubic lattice (b) at n = 1 for t/J = 
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Appendix 

The following high tcnipcraturc expansions up to order have been obtained from the general results. The series 
are expressed in terms oi K = j3J/A. The parameters are X = t/J, z =coordination number, p4 is weak embedding 
constant of square cluster, which has value 0, 1, 3, and 12 for ID, square lattice, simple cubic lattice, and BCC lattice, 
respectively, and ^2 = 2(4p4 + 2; + 2;^)/3. 

Fugacity 



c 



n(l — n) 



2-n 2-n 



--(11 - 27n + ^n^)K'^ - 4^(4 - 2,n)\^K'^ 
8 

-32^(1 - 6n + Zr?)\^K'^ - 32^2(1 + n - n^)\''K'^\ + 0{K^) 



(25) 



Internal energy: 



ujJ = -^n(2-n) (yi^Mz)?)K^VlK'^ -{A-\%n^^r?)K^ -Mzy^K^ 

-256^(2 - 6n + 2,n^)\^K^ - 2h<ozin{2 - n)\^K^ - 5(4 - 6n + Zr?)K'^ - mzX^K" + OiK'') (26) 



Specific heat: 



C^jkn = n{2-n)K'^\{^+8zX^) + 3K -l{'l-18n + 9n'^)K'^ -24zX^K^ 
L 2 8 

-96^(2 - 6n + 3n^)A^/s:^ - 9602^(2 - n)X'^K^] + 0{K^) 



(27) 



Magnetic susceptibility: 

/J-^Xs = ^ + ln{2-n)\l-2K-j{8-6n + 3n^)K'^-4:zn{2-n)X'^K^ 
4 8 L 4 



1 



(2 - 6n + 3n^)i4'^ + 8^n(2 - n)X^K^ 
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+ ^(80 - 180n + 198n^ - lOSn^ + 27n*)iV'^ + ^z{12 + Wn - 41n^ + 36n^ - 9n'^)X'^K^ 
+16zn{2 - n)(7 - 18n + 9n^)X^K^ - 32z2n{2 - n)(l - 4n + 2n^)X'^K^] + 0{K^) 

Charge susceptibility: 



\n{2 - n) - n\2 - nf \{l + 2z\^)K^ + -K^" 

- — (10 - 18n + %n^)K^ - z(l + 6n - -in^^X^K^ 
16 

-8^(7 - 18n + 9n^)\''K^ + 16^2(1 - 4n + 2n^)\^K''\ + 0{K^) 



(28) 



(29) 



ACKNOWLEDGMENTS 



This work forms part of a research project supported by a grant from the Austrahan Research Council. The 
computations were performed on an AlphaServer SC computer. We are grateful for the computing resources provided 
by the Australian Partnership for Advanced Computing (APAC) National Facility. 



12 



* j.oitmaa@unsw.edu.au 



t 



w.zheng@unsw.edu.au; tittp : //www.phys .unsw . edu. au/^zwh 



^ W.H. Zheng and J. Oitmaa, Phys. Rev. B (preceding paper). 

^ H. Roder, R.R.P. Singh and J. Zang, Phys. Rev. B 56, 5084(1997). 

^ N. Shibata, B. Aininon, M. Troyer, M. Sigrist and K. Ueda, J. Phys. Soc. Japan 67, 1086(1998); N. Shibata and H. Tsunetsugu, 

ibid. 68, 744(1999). 
" K. Haule, J. Bonca, and P. Prelovsek, Phys. Rev. B61, 2482(2000). 
^ N. Elstner, and R.R.P. Singh, Phys. Rev. B 57, 7740(1998). 

J.A. Henderson, J. Oitmaa and M.C.B. Ashley, Phys. Rev. B 46, 6328(1992). 

Z.P. Shi, R.R.P. Singh, M.P. Gelfand and Z. Wang, Phys. Rev. B 51, 15630(1995). 
* H. Tsunetsugu, M. Sigrist & K. Ueda, Rev. Mod. Phys. 69, 809(1997). 

^ A.J. Guttmann, in "Phase Transitions and Critical Phenomena", Vol. 13 ed. C. Domb and J. Lebowitz (New York, Academic, 
1989). 



13 



TABLE I: Series coefficients for the internal energy u, specific heat C„, magnetic susceptibihty Xs and charge magnetic 
susceptibihty Xc at /3J = 1 and electron densities n = 1, 0.5 for the linear chain, the square lattice and the simple cubic lattice. 
Nonzero coefficients {t/jy up to order r = 8 are listed. 



r 


linear chain 
n = 1 n = 0.5 


square lattice 
n = 1 n = 0.5 


simple cubic lattice 
n = 1 n = 0.5 



2 

4 
6 
8 


-1.18727275x10"^ -8.78192521 x 10"^ 
-9. 19083667X 10" ^ -6.90488388 x 10" ^ 

2.02166868x10"^ 1.94156599x10"^ 
-5.81316838x10"^ -6.82352222x10"^ 

1.71426415x10"^ 2.25662666x10"^ 


u/J 

-1.18727275x10"^ -8.78192521 xlO"^ 
-1.83816733 -1.38097678 
1.23945460 8.51133343x10"^ 

-1.24540759 -8.52057857x10"^ 
1.38876432 9.81227389x10"^ 


-1.18727275x10"^ -8.78192521 xlO"^ 
-2.75725100 -2.07146516 

3.11186319 1.97093023 
-5.90066730 -2.95078662 

1.32149047x10^ 5.06962222 



2 

4 
6 
8 


1.44045799x10"^ 1.03910138x10"^ 
7.42355042x10"^ 5.62969188x10"^ 
-5.05380430 x 10" ^ -4.99348745 x 10" ^ 

2.39575311x10"^ 2.94484197x10"^ 
-9.61839414x10"^ -1.34812280x10"^ 


Cv/kB 

1.44045799x10"^ 1.03910138x10"^ 
1.48471008 1.12593838 
-3.16863938 -2.18170564 

5.37297355 3.70460027 
-8.40726263 -6.02787617 


1.44045799x10"^ 1.03910138x10"^ 
2.22706513 1.68890757 
-7.98977684 -5.04707069 

2.57521720x10^ 1.27479126x10^ 
-8.15089318x10^ -3.07784724x10^ 



2 
4 
6 

8 


3.02552920x10"^ 2.88871928x10"^ 
-2.69346503x10"^ -1.34119052x10"^ 

1.33858987x10"^ 7.41933727x10"^ 
-6.15230377X 10"^ -3.94194796x 10"^ 
2.69069927x10"^ 1.92550252x10"^ 


3.02552920x10"^ 2.88871928x10"^ 
-5.38693006 x 10"^ -2.68238104x 10"^ 

8.14693283x10"^ 2.59321663x10"^ 
-1.25846811x10"^ -2.89978905x10"^ 
1.94026824x10"^ 3.56078818x10"^ 


3.02552920x10"^ 2.88871928 xlO"i 
-8.08039510X 10"^ -4.02357157x 10"^ 

2.04250289x10"^ 5.55384870x10"^ 
-5.88972872x10"^ -5.65349137x 10"^ 
1.78779922x10" -3.46409964x10"^ 



2 
4 
6 
8 


4.73182256x10"^ 3.60327615x10"^ 
-2.14752162x10"^ -1.24947446x10"^ 

9.57436401x10"^ 6.51283684x10"^ 
-3.93249095X 10"^ -3.25845479 x 10"^ 

1.51554530x10"^ 1.49996060x10"^ 


4.73182256x10"^ 3.60327615x10"^ 
-4.29504323 x 10"^ -2.49894893 x 10" ^ 

5.88984198x10"^ 2.45225449x10"^ 
-8.44234559X 10"^ -2.81650591 x 10"^ 
1.22035314 3.56665468x10"^ 


4.73182256x10"^ 3.60327615x10"^ 
-6.44256485 X 10"^ -3.74842339x 10"^ 

1.47972167 5.40291243x10"^ 
-4.00301487 -7.06439879x10"^ 

1.16033653x10^ 5.01750366x10"^ 



